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ABSTRACT 



In the interest of developing more thorough methods 
of analyzing ship structures, an entire bottom section 
extending from side to side and from one bulkhead to 
another is analyzed as an orthotropic shell of variable 
ciirvature . The longitudinally framed destroyer type ship 
Is studied under hydrostatic loading and longitudinal 
bending stress. Specific consideration is given to the 
orthotropy of the structure and its curvature. 

As a result of the combined analysis of the shell 
and its associated stiffener system, .partial differential 
equations are developed which define the behavior of 
the stiffened shell. The Inherent complexity of the 
developed equations precludes an analytic solution. An 
approximate solution could be obtained, however, by the 
use of difference equations and digital machine computatlcan. 

An application of membrane theory to an equivalent 
orthotropic shell provides a means of determining direct 
stresses and extensional rigidities. A comparative 
study of two existing destroyer types using developed 
theory reveals that membrane theory with corrective 
iterations is of value in prelin4nary analysis. The 
sensitivity of the variation of Hull ciirvature from the 
analytic curvatixre is also demonstrated. The use of 
the midship section coefficient as a design parameter 
for the structural designer should be encouraged and 
exploited. If the midship section area is kept constant 
the curvature can be varied to the advantage of the 
structure with no effect on speed. The relatively high 
bending stiffness in the glrthwise direction of current 
designs merits a reevaluation of transverse strength 
requirements. 

Structural design procedures employing ortho tropic 
rigidity coefficients in combination with either the 
solution to the above differential equations or a modi- 
fied membrane theory, have been outlined by the author. 
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NOTATION 



transverse web frame spacing (in.) 
longitudinal stringer spacing (in.) 

i_ = ( )' 




-Sja° 

external load (Ib/ft^) 

hydrostatic load (Ib/ft^) 

radius of curvature (ft) 

shell thickness (in.) 

equivalent shell thickness (in.) 

axial displacement 

girthwise displacement 

radial displacement 

longitudinal coordinate (ft.) 

distance from shell middle surface 

cross sectional area of a stiffening member (in.^) 

beam of ship at design waterline (ft.) 

dist. of centroid of stiffening member from 
shell middle surface (in.) 

extensional rigidity (Ib/in) 



V 



E = 

G 

H » 

I 

J 

K 

L 

M - 

N = 

Q 





- 

^ - 
f - 



modulus of elasticity (psi) 
shear modulus (psi) 

draft of ship at design waterline (ft») 

moment of inertia of stiffening member about 
its N. A. (in^) 

polar moment of inertia (in^) 

bending rigidity coefficient (lb - in) 

length of shell between bulkheads (ft,) 

length of ship (ft.) 

moment resultant (Ib-ft/ft) 

force resultant (Ib/ft) 

transverse shear force resultant (Ib/ft) 

rigidity moment coefficient (lb“in) 

transverse vertical edge loading on shell (Ib/ft) 

longitudinal vertical edge loading on shell (Ib/ft) 

shear strain 

llneat strain 

Poissons ratio 

stress intensity (psi) 

elastic yield stress (psi) 

longitudinal hull girder bending stress (psi) 
shear stress (psi) 

girthwlse or circumferential coordinate (deg) 
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I. INTRODUCTION 



This analysis can be considered as a contribution to 
t|ae continuing effort being made to produce more rational 
apd less empirical ship structural design procedures. The 
analysis of a structure whose loading and response character- 
istics are not readily defineable does indeed present a 
challange . There exists, however, opportunities for improve- 
ment in the static, elastic analysis of ship structures. 

Such improvements could contribute towards a better determina- 
tion of adequate but' not excessive scantlings arranged in the 
most efficient manner. Scantlings are normally selected on 

basis of a series of approximate analyses. While the effect of 
continuity is often considered, the effects of adjoining members 

and shell curvature are neglected. Checks are made on the 
adequacy of stiffeners and stiffener-plating combinations on 
a panel basis. This artificial method of analysis, while 
adequate and sufficiently approximate in the past, cannot be 
cpnsidered as satisfactory for the future. The knowledge of 
ship seaway loadings is increasing. Advances are being made 
in dynamic loading theory. If these accomplishments are eventually 
to be used to greatest advantage in design, they must be used 
in conjunction with a realistic static analysis. 

The structural analysis of a plate or shell stiffener system 
can be approache'd in general by either of two basically diff- 
erent view points. The first is that of a grillage beam analysis 
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wherein the plating is taken to act as additional flanges to 
the stiffeners in an orthogonal network. Essentially the 
gplllage beam concept as originally developed by Vedeler ^i;* 
rplys on equating deflections at stiffener Intersections and 
moment distribution. The solutions may be obtained by ortho- 
dox solutions to differential equations, trigonometric series, 
energy methods or relaxation and numerical methods suitable 
fpr digital machine computation. Variations in the grillage 
bpam technique and its solutions have also been developed (2). 
One need only to consult contemporary technical journals re- 
lating to aeronautical and civil engineering to find a prolific 
assortment oftechniques and methods for grillage beam solutions 
most of which lead ultimately to machine computation (e.g. 
stiffness matrices). 

The second basic approach is to utilize plate or shell 
t|ieory for the given structure in such a manner that the 
stiffener system effect is incorporated into the medi\am. This 
ip accomplished by substituting for the actual shell with 
stiffeners (frames and longitudinals) an equivalent shell with- 
out stiffeners which exhibits different elastic properties in 
the longitudinal and circumferential directions (orthotropic 
bphavior) . The manner in which the equivalent orthotropic 
shell is determined will be discussed shortly. Let it suffice 
at this point to state that the result of such a substitution 
will be to obtain a shell structure whose "general behavior 

5 Numbers in brackets after an author's name indicate 
references listed in Appendix E of this paper. 
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tqider loading is equivalent to that of the stiffened shell. 

The orthotropic behavior ofthe equivalent shell may be. called 
constructional since it does not result from a property of the 
material but from the manner in which it is distributed. It is 
interesting to note at this point that this orthotropic approach 
is diametrically opposite to the grillage beam method (i,e.,the 
grillage beam assumes plating to act as additional flange 
material for stiffeners whereas in the orthotropic shell, the 
stiffener influence is distributed to the shell or plating) . 

It is felt that the orthotropic shell concept has svifficient 
merit to be the general' approach used in this analysis. Physi- 
cally it appears to be more meaningful in the response and 
behavior of the stiffened structure. While the task of analysis 

of such a structure is not simple, it is at least amenable to a 
mpithematical investigation. Orthotropic shell theory consistent 
with the structure being studied will be utilized. 

Additional methods of shell or plate stiffener systems 
Include the concept that the stiffener system acts as an elastic 
support for the loading and as such may be treated by an applica- 
tion of the beam on an elastic foundation method (3). Another 
interesting technique is to disseminate the effect of one set 
of stiffeners (frames) into the plating while other (longitu- 
dinals) set is subjected to a harmonic analysis producing influence 
coefficients for deflection or bending moment for each 
longitudinal (Ij.) . 

This analysis will investigate the structural behavior of 
the bottom section of a ship from side to side and from bulkhead 
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to bulkhead. The ship type considered will be a longitudinally 
framed naval warship in the destroyer category. JiSraphasis will 
be given to shell curvature, continuity of the shell between 
bulkheads, and the stiffener characteristics. Normal hydro- 
static loading combined with the overall longitudinal bending 
stress will be considered. Shade (5) develox)ed a bending, 
c|aeory for ships bottom structure utilizing orthotropic plate 
theory. However, a flat bottom was assumed and no considera- 
tion was given to hull curvature. 
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II. ANALYSIS 



Preliminary 

The portion of the ships structure to be analyzed is 
sjiown in Figure 2. u & ^ #<> , 
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Fig. 2 Identification of Section Analyzed 

The boundary conditions for this section are discussed 
ijti Appendix C . In brief they are approximated as follows : 

The degree of fixity can be assumed to be 100% clamped at the 
jimction of the shell and bulkhead and simply supported at the 
jtmction of the curved shell to the vertical or flat portion 

ojf the ships side. While these assumed conditions are not in 
fpality present, a determination of the elastic degree of 
fjixlty is outside the realm of this analysis. The postulated 
cpnditions will be sufficiently accurate for approximation 
purposes . 

The inclusion of the huil curvature in the analysis can be 
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effected by the selection of an analytical expression for the 
radius of curvation (r), as a fimction of girth angle {0) . 

T)ie selection of the expression for the radius of curvature 
must be uncomplicated enough so that it will not impede the 
mp-thematical operations involved in the analysis. At the 
same time it must sufficiently describe the hull form for 
any given B/H ratio. A consideration of pertinent hull forms 
and analytical expressions is made in Appendix A . The result- 
ing analytical expression satisfying the above requirements 
hjas been found to be r = 2Hcos0; the radius of curvature which 
describes a cycloid of depth H and span ^ H, One of the sig- 
nificant contentions of this thesis is that a ship type such 
as a destroyer has a hull form which cannot be treated with 
f^at bottom or flat plate assumptions. Of course if a restric- 
tion is made to only a small panel of the structure, a flat 
surface is more closely approached. This type of panel analysis, 
however, is a rather limited view of the structures behavior 

and disallows for the continuity of deformation. Moreover, not 
only is significant curvature present but it is variable around 

the circixraf erence of the hull. An assumption of constant cur- 
vature or a circular cylindrical shell would greatly facilitate 
the mathematical analysis but would not yield results which 
are representative of the structure. It is felt that the 
variable curvature selected represents sufficiently the type 
of hull forms being considered (see Pigs . 19"23) . 
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Isotropic Shell Analysis 

To begin the analysis the elastic laws for the variable 
ciirvature isotropic shell (ships hull without stiffener system) 
will be developed. The shell element and associated stress 
system are indicated in Pig. 3. 

Equilibriijm of forces provide the following relations: 



X dir t 

'$> (N^) dxrd^ d 0 dx f P 2 j-rd 0 dx c =- 0 






-1 







Px 



0 



0 dir : 

9 (N 0 ) d 0 dx 3 (N^^) dxrd^ - Q 0 dx^ - 

-d <ji ^ 




( 1 ) 



~ 3 {(^) d0dxd0 p^rd0dx - 0 

' ( 2 ) 

i 5:^ “ i Qjy -1- P0 -^0 

^ C)0 ^ 



_ 9 _ 



r dir: 



- 2N0dxdO( 



^ (Qy ) dxrd 0 - o) ( ) d 0 dx -f- i^pdxrd 0 — 0 



or 

1 N0 9 Qx ^ 1 ^Qfl - Pp = 0 
r ^x r c>0 



( 3 ) 



Equllibrlvun of moments give — 
about axis ; 

- dxrd 0 - ^ (M ^)d 0 dx 4 * Q(Tjdxrd 0 - 0 

c>x 3 



or 



o>x r o) 0 



(k) 



about p 0 



axis : 



a 




(M^)dxrd0^ a (M0^)d0dx 

W~ 



Q,x‘ixrd0 = 0 



or 



1 3 - Ox = 0 

X r dp 



( 5 ) 
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About axis: 



“N wrd^dx + NM^dxrd^f - 2 %^dxd0 - 0 
r ^2 



or 

V ■ ”9<x+ i «0:. = 0 

r 



( 6 ) 



Elimination of and by use of eq. (ij.) and -(5) gives the 
'following equations of equilibrium . 

i-% + Px = 0 (7a) 

^ X r 



1 - 1 - 1 = 0 



(7b) 
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•Pt- = 0 



(7c) 



^X0 “ ^0x + i ^x - ® ' ' 

r 

At this point there are ij. equations and 8 unknowns. A study 
of the shell deformation will provide additional relationships. 
In the analysis of deformations the following assximptions 

are made : 

(1) All points lying on one normal to the middle sur- 
face before deformation do the same after deformation (i.e., neg 
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lect deformation due to transverse shears). 

(2) For all kinematic relations, the dist. Z of a 

ppint from the middle surface may be considered as vinaffected 
by the deformation of the shell, but that for all considerations 
of the stress system, the stress in the t direction may be 

considered negligible compared to and 0. 

(3) All displacements are small, i.e,, negligible 
compared with radii of curvature of the middle surface and 
their first derivatives are negligible compared to unity 




Fig, k Displacement of Point A 

U, V, W refer to displacements in middle surface of shell 
in directions x, 0, and z respectively. 

Ua, V^, refer to displacements of an arbitrary point A at 
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a point Z from the middle siirf ace . 
from Fig. 4 : 

U - Z 9_W 

3 X 

Va - V / rj_z\ - Z o) W 

' r / r ^ jZJ . 

Wa - W (assiomption 2) 



The strains at point A can be expressed by the following 
relations from Pliigge (6). 

^ x' c) X 

4 i-(b 

Y . ^ ^ 

x0 ^ X ^ 0 



The strains may now be expressed in terms of middle 
surface displacement by using point A displacement relations, 
-z ^^w 



c ^ 

^ " 'iv 
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A simplification can be made at this point by considering 

the shell thickness to be very. small in relation to the. radius 

at any point. In otherwords z is negligible compared to 1. 

r 

or 

(l + z/r) •=. 1 

_1 _ 1 

r f z r 



The strain relations now become 




z \ 

d0^xj 



Carrying out the indicated operations and remembering that 
r is a fxmction of 0 (r = 2HcoslZ^) , we obtain the following 
relations when the notations ^ _ ( ) * and ^ ( ) are 

employed. 



^ - u« _ z w" 

X 



(8a) 






P 



0 




^ (w tan ^ t w) w 

r2 r 



(8b) 




V' 



u 

r 



2 ^ w ' 
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(8c) 
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The stress resultants of Pig. 3 > are by definition 




N 









-t/2 

t/2 

= + 2 /r) dz 

-t/2 

t/2 




(T^(H- z/r) zdz 



M 



xq! ■ " 



-t/2 



t/2 



0 



(1 + z/r) zdz 




dz 



t/2 

^0x r 

-t/2 

(9) 

t/2 

M(2J r - ^ (T 02dz 

-t/2 



t/2 

^x " - y 0x^^^ 

-t/2 



The terms (1 -f- z/r) goes to 1 however by the previous 
assumption of negligible z/r ratio. 

These stress resultants can now be found with the aid 
of the strain relations ( ) . 



- 1 ^- 



since 



(Tx ^-^2^ 

1-S)2 

(FV — 4-V^x) 

1 -v^ 




_E 

2(1+V) 





After substitution of ( 1(3 ) and ( 8 ) 
following stress resultant relations are obtained where 



D 



Et 

1 -D 



2 



is the extensional rigidity and 



K_ Et- 



12 ( 1 -^^) 



TT is the flexural rigidity or bending 
\ 



stiffness . 



N _D (ru* +-Vv + -v>w) 
r 



N . ^ D (v + w + S> ru‘) 

Y> - ^ 

^x0 ' N^x - ^ ( 1 ■ V ) ^ rv«) 

2 r 



(110 
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(11 con't) 






K (r^ w” "V w 4 ”0 tan (f) w) 



M, 






M 



x(Zi 



K (w + tan ^ w w”) 

p2 



K (1:::^) w' 

r 



These relations represent the stress resultants of 
an Isotropic cylindrical shell with variable curvatiire 
r = 2Hcos|^ . (Note that the shell thickness has not been 

fixed yet but can still be a function of ^). These relations 
exist in the unstiffened ships hull. 

We consider now th e gridwork of frames and longitudinal 




Fig, 5 Element of a Gridwork Shell 
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strain at centroid 



**2 = 



where k 



«! - 



then N£ 






^1 


Using the strain relations of 


“0 = !^ 


V - C 0 /w + 


1^1 


r r^ 


similarly 






X ^ 


b 2 


^2 


1 

C\J 

S 

11 








- ^0C0 e 0 


b^r^ 


^1 




- EA0C0 (v/r 


b^r^ 


^1 


similarly for M^, 


Mx ^ Ml - 


N^Cx < 



w) w 

A- 



where w _ change in cixrvature 
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bir" 



2w” 






b. 



- 18 - 






M' 



12 



X _ 



M21 



twist of a stringer =. 0i _ H12 t>i w*^ 



GJ, 



®1^1 _ %2^1 ^X0t>ib2 b^w' 



^(j) 



GJ^ GJx 



GJ^ w 



b2 r 



similarly, _ GJ^ w' 



npte that M0^ Jy( I02 

"^x ^1 



Svimmarizing the gridwork stress resultants of a variable 
curvatiire cylindrical shell form (r = 2Hcos^ } . — 



b2 

V _ ^ 

Mx _ Elxw 

b2 



(u> 




It 



CxW *' ) 



00 (w -f tan)!^) 4. w/r 



r 



EA^CxU’ 

b2 



EAxCx^ w” 

-V- 

b2 



(12) 
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(12 con't) 





r 



NOTE; Nj^gj and N^^c could also be determined but it will be 
seen shortly that they will not be necessary. 

These relations exist in the ships hull stiffners. 
Superposing the stress resultants of the isotropic shell 
with those of the gridwork gives final stress resultants of 
the stiffened variable curvature shell (or the stiffened ships 
hull) , 




r 



^(j) D^(v + w) 4- D^u> - (w + tan0w) 



( 13 ) 



Mx = Kxw" ^ Ko (w 4 tan0w) - S^u* 




w"_|_ tan0wj - (v -I- w) 
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(13 con't) 



N, 



x0 _ + rv») 



N, 



'0x ^ Dx0 4- K0X w 



"V 



X ^x 



w 



M 



0 _ ^X0 



• / 



Where the orthotropic rigidity coeffecients are defined as 



^X . Et ^ EA^ 
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1--P ^ bp 



- Et EA(^ 



1-T? 



4- 



D Et 

S) - 
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1 ^ 
Et^ 

12(1 -i;2) 



4- 



Edx4AxC/) , 



^x0 = 



_ Et 



2(lfV ) 



^ _ 
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12(1 -v^) 
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E(l0 AmC^v ) 



K Et^ 

'V •=• 



V > _ ^x^x * E0 _ EA0C0 



12(1^) 



b£ 



Ex0 - Et3 



GJ, 



+ 



^7 



12(HV ) 



Et- 



12(1 + V ) 
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GJ^ 

bn 



(li;) 
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Note that the and N0^ represent shear carried by 

t}ae shell alone. The grldwork can only carry this shear 
by a bending of its stiffeners, a bending which becomes 
impossible when the shell is joined to the stiffeners. 

Introducing, now, these new stress resultants for 
the stiffened ships hull, into the equations of equilibrium 
(eqs. 7), one obtains the following differential equations. 




r 



^ (V -^w) D^u' 



S0(w tan^w) 




D^(v 4 -w) , u' - ^(w +tan 0 w) -1 - K^w' ' Kv (w" + tani^w") 

r2 r 

- S^u"|^ 1 ^ w'"j 

r d 0 \ ^ / 






r ^ ^ / r 
V 



r 




tan0wj 






(v + w) 



- P, 



0 
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After carrying out the indicated operations, rearranging 
terms and using partial operator notation as follows. 



d, , ^ ( ) ' 


^^2 ^ 


_ . ( ) 






3 j 2 T 




and d]^^ _ o>^ > 


d2^ ^ 














the differential equations assvime 


the form. 






^ir^d^ + ^2^^^2 ■*“ G2rtan0d2^ 


u + [03 


r^dfd^ V 


( ISa) 






2 Cgtan 0 8,-j^8.2 


]w . -r 3 p. 







+ C0r tan0 v + " CiQr^d^d2 + Cyr^d, 

-3GQr tan20 d2 - C3^2.*^2 " ^ 2 ^ 



(15b) 



- 2 Cg tan0 d£^ - C-]_q_ tan^ d£^ 4 Gy r^tan0 + Cgr tan^ w--r ^0 
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4 ■* ^*^8^ taxi^P d£ - Cgrd^ v -i- 

^12r^di^ -\- + S*^2^ ■<• C^^r^tan0 ^^^^2 

+ Cj^^tanji d2^ - 2,C0rd2^ 4- 8C<^tan^0 <^2^ + 2G^d2^ (15c) 

4 2C]_^d2^ 4 7C2_2_ tan^0 d 2 ^ - 4Ggrtan0 d 2 4 9Cxitan0 d 2 
4- l5Cxi tan^0 d2 4 Cyr^ - 3Cgr tan^0 - <-0^ w •=. 

The constant part of the coefficients are identified as 



°1- 


Dx 


°6 = ^0 


°2 = 




09=1^ 


O3 = 


^-0+ 1^X0 


®10 = ■*■ 


11 

0 


I'.) 


Oil. 


C5 


Sx 


C12 - 


“6 = 


K0X 


^13 = ^ X 0 ■*- 


° 7 ' 




= 5 *^ + 



hi 



Equations (15» a, b, c) comprise a set of 3 partial differential 
equations in 3 unknowns, (u,v,w). Their solution as a function 
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of radius and rigidity coefficients will completely define 
the structural behavior of a stiffened cylindrical shell 

(of cycloid cross section) for given loading and boiindary 
conditions. The rigidity coefficients were made constant 
as a necessary simplification. It will be seen later, how- 
ever, that many are in actual practice essentially constant 
over considerable portions of the shell section. Additional 
comments concerning these differential equations will be 
delayed luitil the Discussion of Results section. 

Membrane Analysis 

Since the solution -of the developed differential 

equations is not readily available, one must approach the 
shell problem with other and more immediately useful tools 
of analysis. It must be stated at this point, however, that 
the analytical approach by means of the differential equations 

of the problem is the only means by which a complete and 
relatively exact solution may be obtained. It is unfortunate 

that their complex nature is currently limiting their use 
by shell designers. The alternate methods employed to 

o)Dtaln solutions include the beam method, iterative calcula- 
tion techniques, membrane theory, and other unique approxima- 
tion approaches. Lundgren (7) describes rather completely 
the various methods of shell analysis currently in use. 

The applicability or the relative desireabillty of the 

y 

alternate methods is determined basically by the shell dimen- 
sions, The ratio of the shell length to the shell span (or 
beam), the L/B ratio, normally will govern the technique. The 



- 25 - 



application of membrane theory (i.e., the consideration of 

direct stresses only; the support of the load by tension, 

cpmpression or shear forces only) is considered good for 

short shells or shells whose L/B ^ 1. Lundgren (7) states 

that the membrane condition of stress is a very good 

approximation for the short shell since the elastic support 

between end sections (or bulkheads) is quite stiff and 

loading is balanced by N N, and N ^ stresses. The 

^ jO yi(p 

ASCE Manual of Engineering Practice No. 31 recommends a 

superposition of membrane stresses and those created by 
boundary conditions. However, on the matter of boundary 

conditions, the manual also states that for a L/B < .83 

the transverse (Bulkhead) end effects damp rapidly and the 
effects of the longitudinal edges are felt only in a small 
area. 

The previous discussion indicates that a membrane 
analysis of a ships bottom section between bulkheads is a 

reasonable approach when one realizes that the L/B range 
for the ships or ship type under consideration is as indica- 
ted in Table 1. 



DE 1033 


DDG 952 


DLG 9 


DLG(N) 


CG(N) 


CA 134 


0.775 


0.85 


0.81 


0.965 


0.56 


0.56 



L = length between bulkheads ( \m ^ ve<^ic.^r) 
B = beam at waterline 



TABLE 1 Length to Span Ratio of Ships Considered 
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The membrane stresses (Ib/ft) in the shell element can 
be obtained by setting the shell moments equal to zero in 
the equations of equilibrium (eqs» 7). 

“fj - Pi.=^ 

”c) ^ “Pp " ^ ^^0 

^ X r 'd0 

^ -Px - 1 

^ X V gf 



integrating the latter two of these. 





0+ i 

r bp / 



dx + 1 ( 0 ) 



N 



X 





dx -I- f 2 (^^) 



for =® 0 and and p^^ Independent of x. (These are correct 
conditions for ships hull loading. The longitudinal stress 
resulting from the hull girder action and determined from 
the longitudinal strength calculation will be accounted for 
shortly but will not be considered as a p^ loading at this 
point) . 

^x^ = " (P0 ^ £ (d ^0^ X 4- f]_(0) — “xP^0j -)-• f-j^(0) 

^ ^ UiVwtfi. + Jj. 

Nx = dF(0) - X dfi(0) ^ f2(0) 

2 r d0 r d0 
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P(0) is determined by loading and shell shape 
(oxxrvature) but i*i(0) and f 2 ( 0 ) will be determined by 
bp\mdary conditions at the shell ends (bulkheads). 

Consider first the case of simple support end 
(bulkhead) conditions. 







y 


X' 


^l<eel 

ir 


At X - t L/2 , 


o 

II 




- '■k- 


- - 




then. 




Pig. 6 


Shell 


on Simple Supports 



l 2 dP(0) - L df]_(0) f£(0) 

^ 8r d0 2r ^ d0 






8r 



dP(0) L 
~d0 



df3^.(,0) ^f£(0) 



2r d0 



which gives, 

fl(0) 0 ^2(0)= - dF{p) 

8r d0 



arid therefore. 




- xP(0) 



^ -_1 (L^ - ipc^) 

8r 



dP(0) 

"d0 
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Now consider shell end or bulkhead conditions for the 
shell continuous over the supports and assume 100^ clamping. 
As an aid in determining the stress over the supports, 
ope can make use of the relationship given in ASCE Manual 
No. 31: the ratio of the stress (in a shell over continuous 

supports) at the point of support to the stress in simply 
supported shell at the midspan is equal to ratio of the 
moment in the continuous shell at the support to the moment 
in the simply supported shell at the midspan. 



ov- 

_ % 

^ = “2/3 

NxE - - 2/3 N^q 

and N^q was foimd previously 

N at support (X -tL/2) 
yi 

continuous over supports. 

At X = i L/2 

_ q _i£ dF(/6) 

12 r d0 



and, for 

“o = 



a strip of shell, 
-1/12 wl2 



1/8 wL^ 



£ 






Fig. 7 Comparison of Continuous & 
Simple Supports for Shell 
to be - l 2 dF(0) 

8r d0 



is 1 l 2 dF((2() 



I2r d0 



for a shell 



Vx - I - 



Ice®*" 






Fig. 8 Shell Continuous over Supports 



or 



dF((Z() 

d0' 



l 2 dP ( 0 ) 

12r d0 8r 



L dfi( 0 ) ^^ 2 ^ 0 ) 
2r dj?^ 



dfl(0) 
d 0 

wjiich gives, 

fl i<f>) ~ 0 

f2 (</) - - _l£ dP(0) 

2I4J? d0 



l 2 dP(0) 
12f d0 



8r 



dP( 0 ) L 

-\r 

d0 2r 



then, 




xF( 0 ) 

(L^ - 12 x^) dF( 0 ) 

21+ r d0 



^0 = Pr ^ 

where , ' 

P(0)= (P0^ 1 'ct N^ ) 

r d 0 



( 16 ) 



Equation 16 represent the general form of the membrane 
stresses in the ships bottom section of radus of curvature 

r . 

These stresses will nox'j be developed in a more final 
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form by considering the hydrostatic loading and the 



ction cylinder. 



gig, 9 Loading on Ships Hull 
0^ ^ Y (Vxo + ^) 

_ XL*"" 

Pp * -(q - w cos 0) 
f^r = " ^ C ^0 + ^ 1 cos20 + w cos ^ 

, = - w sin 0 

-Y r ^0 ■*' H/2(1+ cos20)J 4- w cos 0~j 2H cos p 

_ - Y [ 2H hQ cos 0 + 2H^ cos^ p] 4 2wH cos^0 ^ (7 

and since P(0) •=. P0 1 d ^0 

r d 0 

P(0) = Y j^hQ tan0 4 ^ sin20] - 3w sin0 

and are then found by applying equation (16) . 

^x0 = - X 1^^ ^^0 fcan0 + 3/2H sin20) - 3w sin0] (17b) 

N L^- 12x2 .y(ho sec^0_4 3H cos20 ) 1 (17c) 

^ il8H“ L 



Let qo := V hQ where hQ - an 
arbitrary design head to account 
for wave forces on freeboard surface 
Of hull. 

z = H/2 (I4cos20) from appendix (.A ) 

w = internal loading of shell 

p 

due to machinery weights (Ib/ft^). 
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Equation (17) are the membrane stresses for the ships 
bottom section. With the aid of these relations, the necessary- 
shell thickness (that is, the equivalent shell thickness for a 
hull with no stiffeners) can be foxmd as it varies arovind the 
girth. From this equivalent thickness, extensional rigidity 
coefficients (D^ , ) may be computed. For the 

i” i“ ~E~ i 



single equivalent thickness derived there results the fact 



that Djj. ^ since they both equal t@ or l.lt^ see appendix D, 
i E 

However, in practice these coefficients are very nearly equal 
as can be seen from Figwe 26. 



Tx 



N 



12 t. 









12 t, 



( 18 ) 









12 




Figure 10 Direct Stresses on Equlvalen-|: 
Shell Element 



Where = the hull girder longitudinal bending stress as 
determined from the longitudinal strength calculation. 

Using the maxim-um shear failure theory and in particular 
the Tresca version, 



rr- 

_ I max 






+ r-'p 



or 



(19) 



- (Tx^ - ^ 
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Substituting equation (l8) in equation (19) and 
solving for t^ gives^ 

t^ _ - B \Jb^ - i^AC 

2A (20) 

where A ■=- ll|4 ( 

B - 24 ^i(N^ - ) 

C = n/ - 2N^N^ 4n/^ 

It now becomes necessary to make a judgement as to the 
treatment of and ';p , 

y 

d” will be assumed as 33»000 psi and a safety factor of 
1,65 will reduce its value to 20,000 psi, 

(J"^ will be approximated by Tobins relation that 
(T- where Lg = length of ship in ft, 

r Ton/ in^ 

If it is desired to have ^ as a ftmction of ship depth 

(or girth), with - d~-i Q-t the keel (0 ^ 0) and 

1 1 max 

_Ti max at the waterline or N.A. (0 ^2) (to allow 

for the ship being in the heeled condition), 

then Lg' (5 -h 3 cos20) (21) 
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The membrane analysis has now been sufficiently developed 
to make a test of its applicability to specific examples. For 
a given ship (or design), the equivalent shell thickness 
variation can be found by means of equation (20) which requires 
a determination of membrane stresses equation (1?) and hull 
girder bending stress estimate equation (21). The extensional 
rigidity coefficients are then readily calculated since the 

equivalent shell has no stiffeners and ^ t@ , 

E E 1- 

T)ie calculated values may then be compared to the actual 
values on the ship selected for study. It should be noted 
that in the calculation of the membrane stresses, the longi- 
tudinal section (X — constant) selected for stress determina- 
tion was the section at the bulkhead. This was a logical 
design decision since a close examination of equation (20) 

(or a few comparative calculations) will indicate that analysis 
of this section results in the greatest tg for any given ^ 
value along the shell length. In addition, the value of w or 
internal loading was assumed equal to zero. Not only will 
this be the more severe loading condition, but it is also the 
one which current design (series of flat panel analyses) 
employs. This assumption will assist in making the comparison 
more valid. The retention of this quantity in the general 
theory, however, is recommended. Additional comments on 
this subject will be made in the discussion section. 



III. RESULTS 



The analysis of the orthotropic shell by superposition 
of the isotropic shell and a stiffener grldwork (with all 
effects of bending and torsional moments in the shell con- 
sidered) resulted in a set of 3 partial differential 
equations in 3 unknowns (u, v, w). They are equations (15) 
op page ZZ . Due to their lengthy nature they will not be 
repeated here. The coefficients in the partial differential 
equations are variable in that they contain the radius of 
curvature in powers up to the fourth. Equation (l^c) is of 
the fourth order while l5a & b are third order. Equation (l^a) 
could be made homogeneous by setting p^ = 0 but equations 
l^bjC must remain non- homogeneous , All three partial 
d^fferentiai equations are linear. 

The membrane analysis produced more tangible results 
in the form of expressions for the membrane stresses as a 
function of: girth angle (p), position along shell length (X), 

and ship characteristics such as ship length (Lg), draft (H), 
compartment internal loading (w), distance between bulkheads 



a A bottom section constructed with no stiffeners but an 
equivalent thickness (tg) to produce the same extensional 
rigidity as one with stiffeners. 



(L), and a design head on freeboard plating (ho). 



Membrane stresses in ships bottom section .'^- 
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Figure 11 Illustrates the membrane stress distribution 
along the length of the shell (a section const). 

Figure \X represents the circxamf erential distribution 
of membrane stresses' in i)E 1033 and DD 931 hulls . (a ^section x const). 
The derived expression for the shell equivalent thick- 
ness variation utilizes the membrane stresses along with 
the hull girder bending stress ( (Ti) and a yield or 
allowable design stress ( (Ty) • 

tg ^ - B ^ \)b^ - lpA.C Inches 

— 

where 

A = 141(. ( . (P^2) ^ 

. B = 24 (Ti (N^ - Njj) membrane stresses = Ib/ft 



by 



The extensional rigidity 



^5 

E ” E 1- ^ 



1 ,lte 




tg __ » 38 tg 

2(1+V) 



D 

V - 
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coefficients are then computed 

for 0.3 



Using the DE 1033 and DD 931 hulls, these coefficients 
have been calculated from the developed expressions above and 
compared to those calculated from the actual hull structure. 
The results are plotted in Figuresl4-» ^5, 16, 17, 
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Pigvire 11 Longitudinal Distribution of 



C>i 




37 



G 1 T>\ST^ qT Sri^g^Sgs it^ SH\P^ I-\ul.l, 




{ 









38 » 



Ci- ^ 



T“ 




I 



Figure ll|. Extensional Rigidity Comparison 

for DE 1033 
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Figiire 16 Extenslonal Rigidity Comparison 
for DP 931 
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IV. DISCUSSION OF RESULTS 



The inherent complexity of the developed set of 
partial differential equations precludes a simple solution. 
An analytical solution appears to be a forraidable challange. 
Perhaps the most promising course of action would be to 
apply the technique of difference equations and with the 
aid of a digital computer obtain an approximate solution. 

The complex nature of the equations should not discourage 
efforts to solve them, for in their solution lies the 
complete behavior of the stiffened shell. This last state- 
ment needs qualification in that the behavior does not 
include the effects of concentrated loads, or localized 
stresses. However, as noted and discussed in Appendix B 
the general behavior of the shell should be obtained first 
and then a follow up check made for localized stresses. 

The shell analysis should not proceed solely on a series 
of localized stress analyses throughout the structure. 

The value of the above equations is naturally dependent 
upon the validity of the included assiamptions . Perhaps the 
most questionable one of these is the assumption of con- 
stant orthotropic rigidity constants with variation in 
This simply states that the structure is uniformly con- 
structed of the same size structural elements around the 
section. This, of course, would never be done since the 
loading varies with j/. Figure 26 illustrates how the 
actual coefficients do vary with It can be seen that 



there are regions of the shell where many are constant 
with 0. For those that vary slightly, a mean value might 
be taken. The sharply varying ones are associated with 
the stlffners in the ^ direction. Perhaps a reanalysls 
of the transverse strength requirements will reduce 
these to more smoothly varying values with 

Once the equations are solved, the solutions will be 

in the form of expressions for the deformation as a fxmc- 

tion of X, j() , and the rigidity coefficients at any given 

point ^with, of course, the effect of the injected boimdary 

conditions being incorporated in the functional expression). 

Once these deformations are known, expressions can be 

obtained for the strains, moments, shears and stresses in 

the equivalent orthotropic shell. The rigidity coefficients 

may then be determined at any point by assignment of design 

stresses. These design stresses would, of course, be 

influenced by considerations of hull girder bending stress 

and stability of the structure. After all the rigidity 

coefficients are determined (i.e., the structural member 

sizes and spaclngs) they could be fed back to the original 

expression for deformation along with their variation in j6^ 

A recycling of this process until the coefficients do not 
change appreciably will give the final structural design. 

A follow up check for localized stresses and concentrated 

loadings will complete the process. 

The application of this theory and procedure to a photo- 
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elastic scale model and later a larger steel structural 
model would prove invaluable. A program of testing, 
modifying theory and' reorienting basic approaches would 
most probably result in new and improved concej^ts in 
naval structural design. 

A few comments are in order concerning the hull 
curvature. The reasons for selecting the analtical 
curvature (or radius of curvature) used in this thesis 
are ai» 5 )ly covered in Appendix A. However, for additional 
research in this area, NACA TM 1302, pp, 8 - 16 is 
strongly recommended. The topic of variable curvature is 
discussed and fourier polynomials are developed for the 
radius of curvature expression. Very simple expressions 
result from the analysis and they could be used to approxi- 
mate, very closely, any given hull design. It is quite 
interesting to note that within the scope of preliminary 
hull design, there exists a design variable xmused by the 
structural analyst. This design parameter is the miship 
section coefficient (the ratio of the immersed area of 

the midship section to the product of beam and draft, Ajg ). 

BH 

This coefficient defines the relative fullness at the 
midship section to a circumscribing rectangle. In effect 
it is a shape factor which determine the shell section 
profile. If the stress analyst had some degree of control 
over this coefficient, a shape or section profile could be 
determined which would be the most beneficial for the stress 
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in the hull. The stress magnitude and distribution in 
any shell structure is a function of its shape or curva- 
ture. The encouraging fact is that indeed the analyst' 
can play some part in determining this coefficient to 
the advantage of the structural requirements X'jithout 
affecting the resistance of the ship, D. W. Taylor (8) 
demonstrated this fact in 1908 when he conducted model 
experiments to determine the influence of midship section 
shape upon the resistance of the ship. His tests were 
definitive enough for him to make the following conclusion; 
"for vessels of usual types and of speeds in knots no 
greater than tv;ice the square root of length in feet, the 
naval architect may vary widely the midship section 
fullness with no material beneficial or prejudicial effect 
upon speed." It must be remembered, of course, that the 
above statement applies to the shape and coefficient of 
the midship section for a given area and not to the area 
itself, Taylor varied the midship section coefficient 
from 0.7 to 1.1 using 10 models and holding the midship 
section area constant on all tests. Thus it appears that 
the structural analyst can accept a given design and have 
considerable freedom to determine his optimun shell shape 
(for a given midship area) within the limits of not adverse- 
ly affecting the pin s c coefficient or general 
seaworthiness . 

The use of membrane theory for the evaluation of stress 



I 



in the ships bottom section has considerable merit. 

Although the stiffened ship structure will carry the 
stresses according to the distribution of its members,' 
membrane theory will predict what the stresses are that 
must be carried. 'Ihere are two basic limitations to the 
theory. (1) The normal shears, bending moments and 
twisting moments are neglected, (2) the membrane theory 
becomes invalid at and adjacent to the boundaries. These 
limitations are not severe, however, since the quantities 
of normal shear and bending and twisting moments are 
normally of a very minor magnitude. (It will be shown 
shortly how this can be determined quantitatively) . The 
boundaries can be subjected to a separate discontinuity 
analysis and do not jeopardize the evalution of conditions 
over the major portion of the shell. To prove that the 
membrane theory gives is a good representation of the 
situation, one can proceed in the following manner (1) find 
the displacements associated with the direct stresses pre- 
dicted by membrane theory* (2) use these displacements in 
the general equations which incorporate bending (equations 10) 
(3) compute the direct stresses from these bending theory 
equations and compare them to the membrane direct stresses 
to evaluate the accuracy of the membrane theory, 

PigxiresU, 12, represent the membrane stress distri- 
bution in the longitudinal {(f)- const) and girthwise (X- const) 
directions respectively. It must be kept in mind that 
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these membrane stresses are invalid at or adjacent to 
any of the shell boundaries (bulkheads and longitudinal 
junctions with vertical ships side). A matching of 
reaction forces at these locations is necessary as well 
as the determination of secondary bending to satisfy 
compatibility. 

The most Interesting feature of Figure ]J. is the 
presence of longitudinal compression (- stress) in 
the shell midway between the bulkheads. If the ship 
actually had a flat bottom the stress would obviously 
be a tension. As the bottom curvature increases and 
becomes a shell, the longitudinal stress does Indeed 
become compression. The direct stress whose magnitude 
and direction is the most sensitive to shell shape or 
curvature is the longitudinal direct stress, N . This is 
not a conclusion of this thesis but rather a standard fact 
of shell structure behavior (a fact which is largely 
responsible for the widespread success and popularity 
of concrete shell structures). This is of Interest to the 
naval structural engineer since it is this particular 
stress which is combined with the hull girder bending stress. 
The recognition of this stress should render ship structural 
stability considerations more realistic. 

Figure 13 represents a plot of the longitudinal stresses 
as measured on the destroyer HMS ALBUERA (9) in the hogging 
and sagging conditions. The test section (fr. 62 ^) was 
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located in the midship region and about halfway between 
transverse bulkheads. Considering the observed stresses 
in the hull above the waterline, a linear relationship 
is readily established. When the imderwater portion of 
the hull stresses are examined it can be seen that they 
exhibit a variance from the established linearity (the 
linearity established the upper portion of the hull must 
be the same as for the lower since n^her moment nor 
section inertia has changed in going to the lower portion) . 

This departure is obviously due to the effect of the hydro- 
static loading on the underwater shell or hull. It will 
be noted that this effect augments the negative hogging 
stress and reduces the positive sagging stress. In other 
words the hydrostatic loading on the hull produces a 
longitudinal compression This agrees with the theoreti- 

cal plot of N in Figure 11. The longitudinal compression 
effect imder hydrostatic loading would be no surprise to the 
civil engineer who works with shells. It would, however, be 
unexpected by a naval structures analyst who limits himself 
to the study of small panels and predicts only a flat plate 
tension effect due to hydrostatic loading. While it is quite 
true that each panel as such is subjected to a tension pro- 
ducing load, this is only a minor superposition on the more 
general shell behavior which is a compression. These local- 
ized stresses can be checked in the later stages of design. 

Figure 12 illustrates the distribution of direct stresses 
around the girth (X= const) for two selected destroyer types. 
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The section portrayed is at the left transverse bulkhead of 
the basic shell (X= -L/2). This particular section was 
selected for purposes of obtaining a design shell thickness 
(and its variation with 0) since its thickness requirement 
is greater than any other intermediate section for any 
given ^ value. This was an unfortimate choice, however, as 
a forthcoming discussion will indicate. The direct stresses 
of Figure 12 were used in calculating the equivalent shell 
thickness and hence extensional rigidity constants (tables ij. , 
6, ). These coefficients are compared with those actually 

existing in the ships as a result of their construction. As 
was previously stated, the membrane theory will give reason- 
ably good results everywhere except the boundaries. This 
unf ortimately vjas the necessary region (trans. BHD) to 
analyze since it produced the controlling equivalent shell 
thickness. The evaluation of the direct stresses could have 
been and should have been improved by using the procedure 
indicated on page 14.8 . Since the extensional rigidity coeffi- 
cients were calculated from tg which in turn was calculated 
from the direct stresses; additional corrective action on 
the direct stresses produced would result in a different and 
better correlated set of extensional coefficients. It is 
surprising that at this worst possible location for membrane 
theory, the degree of agreement between constructed and cal- 
culated is still reasonably good. A corrective cycle on the 
obtained membrane stresses would result in a general increase 
of varying degree for the rigidity coefficients at all values 
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of 0 (since all 0 values are on the boundary). The differ- 
ence in relative variation of calculated from constructed 
values in the two ships selected is due to several factors 
which characterize each ship. It is believed, however, 
that one of the major causes is the fact that shell curva- 
ture of DD 931 (considered similar to DD 952) resembles a 
cycloid to a lesser degree than does DE 1033* This sugg- 
ests the importance of selecting the most accurate analytical 
model for hull curvature. Other causes for non similarity 
of the ships plots are possible opposite departures from 
the assumed. conditions in the membrane equations (e,g, 
stress estimate, factor of safety used and degree of fixijty 
at ends), Plg\ires 15>17 > merely represent the comparison 
of the calculated and constructed ship hull thickness since 
these rigidity coefficients are just a product of a constant 
and the thickness. 

The approach used in this thesis was not expected to 
produce any radical changes in structural design. It was 
intended to achieve a more rational design attitude by treat- 
ing a section between bulkheads as a shell and then combining 
these results with the dverall ship bending stress. An 
expected result was a slight reduction in structural members 
as determined by current design. With due regard to the 
previously discussed obstacles and handicaps involved, the 
plots of Figure Ii|.,l5>l6 »17, Indicate an apparent trend in 
this direction. The determination of strucutral members can- 
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not be made on extentlonal coefficients alone. A recycling 
of the membrane solution as noted on p will give magnitude 
of moments in the equivalent shell. The individual structural 
members for the ship may then be determined by a combined 
consideration of the extentlonal rigidities required, the 
moments in the shell and the localized panel stress. 

Unfortunately neither one of the basic approaches used 
in this thesis (l.e., shell bending theory by differential 
equations and membrane theory) has been followed through 
sufficiently to give definitive quantitative conclusions. 

The investigation, however, has been quite profitable in 
that the ships hull has' been analyzed in a manner more 
appropriate to the geometry of the structure. Such an analy- 
sis can and, indeed, did reveal information unobtainable 
by an otherwise less realistic approach. Specific areas for 
additional study and analysis have been clearly indicated. 
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V. CONCLUSIONS 



1 . The resulting partial differential equations of the 
orthotropic shell behavior are complex and even after 
simplification are not amenable to an analytic solution. 

2, A rational structural design procedure can be 
established by using the partial differential equations 
solution and the orthotropic rigidity coefficients (see 
page 4-5" ) . 

3« The analytical expression for the hull radius of 
curvature must be uncomplicated and yet representative 
of the actual hull form being analyzed. 

I 4 .. The midship section coefficient or the midship sec- 
tion shape is a design parameter which can be used within 
reasonable limits to benefit the structural design with 
no prejudicial effects on hifLl resistance or ship speed, 

5 . Subject to its characteristic limitations, membrane 
theory can be used in the preliminary analysis of ship 
structures. Iterative corrective procedures can be 
Used for improvement .of solution and/or estimation of 
existing error, 

6 , A longitudinal compression results from hydrostatic 
loading on a ship hull that possesses curvature. The 
magnitude and distribution of the compression around the 



girth of the ship is sensitive to the hull shape or 
variation of c\irvature. 

7. Useful extensional rigidity coefficients can be 
determined by membrane theory when Iterative corrections 
are included, 

8. Current design methods produce hulls which have a 
girthwise bending stiffness' which is five times as 
large as the longitudinal bending stiffness. 

9. The different and more appropriate approach used 

in this thesis for the analysis of ship structures with 
hull curva^ture has produced; 

(a) A specific knowledge of items within this 
analysis requiring ftirther study and refinement, 

(b) A more discerning attitude towards the 
structiire being designed, 

(c) Suggested design procedures utilizing the 
methods and techniques discussed in this study. 



- 5 ^- 



VI. RECOMMENDATIONS 



1, Computer solutions to characteristic partial 
differential equations of the orthotropic shell structure 
should he devised, programmed and tested for merit in 
relation to cxirrent or other proposed design concepts, 

2, Optimization of ship structural members and their 
arrangement should be on the basis of an orthotropic 
shell and not on a series of panels or current design 
approximations . 

3, A study should be conducted to determine if a 
defineable relationship exists between the longitudi- 
nal and transverse strength requirements of any given 
design. The current ratio of five to one on girthwise 
flexural stiffness to longitudinal stiffness appears 
excessive . 

Ejq^erimental verification of all developed theory 
and techniques should proceed conciorrently to aid in 
providing realistic results which are of a practical 
significance. Experimentation could commence economically 
by photoelastic stress analyses of small scale models. 
Tests on successively larger steel scale models using 
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adequate instrumentation should proceed if early resiilts 
indicate that further testing can be of value, 

5, A short study should be made to determine an analyti- 
cal expression for the radius of curvature of any given 
design. The expression must be uncomplicated or capable 
of being simplified for the mathematical operations 
which must be made on it, (NAGA TM 1302 pp. 8 - 16 is 
recommended for this study) , 

6, Since the mathematical analysis of a circular 
orthotropic shell is quite complete, one possible approach 
is to break the shell into a finite nvimber of segments 
represented by constant radii of curvature (circular), 
make the analysis and join the segments by making them 
compatible at all jxmctions, 

7, A strongly recommended topic for study is the 
variation of midship section coefficient or shape for a 
given midship section area and the determination of the 
optimum shape for structural design. The only major 
constraints will be the maintenance of the designed 
prismatic coefficient and desired seaworthiness character- 
istics. Within these constraints, however, there is 
latitude for increasing structural efficiency. 



8 . 



A modified membrane theory should be developed for 



the preliminary analysis of ship structures. Since the 
solution improves with iteration, it is quite amenable 
to machine computation, 

9. A method should be devised for readily determining 
the magnitude of longitudinal compressive stress 

which is produced by hydrostatic loading on all hulls 
with curvature. This value should then be included 
in any deliberations on structural stability. 

10. A complete analysis of the stability of the 
ortho tropic shell between bulkheads shoiald be combined 
with this thesis to provide the complete ingredients for 
a new concept of ship structural design, 

11. A relieable solution to the partial differential 

equation of the orthotropic shell will require a know- 
ledge of the^ shell boundary conditions. It is apparent 
that a study devoted soley to this subject would be 
of considerable value. 
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A. HULL CURVATURE 



The curvature of the hull in the girthwise or circum- 
ferential direction is one of the most significant factors 
in this investigation. The author feels that the magnitude 
and variation of hull curvature for ships of the destroyer 
type has a significant effect upon the structural behavior 
and therefore should be given full consideration in the 
structural analysis of these hulls. The section of the 
ship considered was in the midship section region. There 
will normally be no longitudinal or fore and aft hull curva- 
ture in this middle portion of the ship. 

The argument could be made that the hull girthwise 
curvature is really negligible if a basic panel (botinded by 
frames and longitudinal stringers) is selected for analysis. 
The weakness of this argumenl^, however, is that the hull’s 
general response is not on a panel basis but on a much larger 
and continuous region (i.e., the portion of the hull structure 

between transverse bulkheads) wherein any flat plate or 
flat bottom assumption would be xAnrealistic . Additional 
comments on the relative size of hull portion analyzed will 
be reserved until Appendix B. 

The task of introducing the hull curvature of a specific 
type of ship, or of a given design within a type, is a diffi- 
cult one. The analytical expression of curvature selected 
must be representative of the design being studied. It must 
also beAreasonably Aoncomplicated function capable of being 
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subjected to the necessary mathematical operations without 
excessively expanding the complexity of the resulting 
analysis. The representation of ships hulls by polynomial 
expressions is not appropriate since the desired expression 
is one giving the radius of cxirvature as a function of the 
girth angle (^). 




The analytical expression finally selected as satisfying 
the above requirements was the equation of the cycloid 
r =r i|.a cos (j) (Figure ) where the depth of the curve 
(draft) is 2a and the span (beam) is 2lTa. ("a” is simply 
the radius of the generating circle). This type of curve 
oj* section profile results in a span to depth ratio (or B/H 
ratio) of 'IT . The B/H ratios of the ship type considered 
is indicated in Table 2. 
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DE 1033 


DD 952 


DLG 9 


DLG(N) 


CG(N) 


CA 134 


B/H 


3.16 


3.13 


3.14 


3.02 


3.25 


3.14 



TABLE 2 . B/H Ratios of Ship Types Considered 



The 'll B/H ratio of the cycloid does not in Itself 
qualify this curve as being a good representation of the 
destroyer hull. Tlie curvature and its variation with 
must correspond sufficiently close to the ship hull shape. 
This is the most significant factor. It must be remembered 
that coincidence of the analytical expression and the ships 
hull is not of importance (i.e., the r, )6 relationship should 
not be misconstrued to be a polar plot). The relative 
'equality of the values of the radii of curvature at given 
0 values is the significant measure of the expression's 
adequacy . 

The parametric equations of the cycloid in 0 are, 

y ■= a(2^ -t- Sln2)6) (22a) 

z = a(l + 2^ (22b) 

The cycloid section can then be fitted to any given design 
by letting a = H/2 and plotting 



y 




obtained from 



an elimination of (j) in equations 22. These plots have been 
made In figures 21 , 22 , 23 , along with the plots of the 
curve produced by fitting an elliptical section to each actual 
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hull shape selected for comparison. 

A study of the comparative plots of hull shape and 
analytical expressions reveals that the elliptical section 
possesses a greater resemblance to the curvature (or radius 
of curvature) of the actual hull than does the cycloid (note 
how the DE 1033 curvature is virtually coincident with the 
elliptical section) , However, the analytical expression 

for the radius of curvature of the ellipse is 
2k2 

r = where a and b represent the major 

a^sin^0 b^cos^0 

and minor axis respectively. This expression for the radius 
would not be suitable for the mathematical operations in- 
volved in the shell theory of the hull. It is also of 
interest to note how the hull form of the GA 134 is not 
appropriate for this analysis since it deviates excessively 
from curvature given by either expression, CG(N) hull 
form would also be considered unsuited for the cycloid approxi- 
mation. 

The obvious disadvantage of the cycloid section is the 
\mfortunate fact that the radius of curvature goes to zero 
at 0 - ^/2 , This is inherent in the nature of the cycloid 
curve since at the extreme points of the span the generating 
curve changes the sign on the radius of curvature and, at 
this point, it passes through a zero point. This is not a 
major disadvantage, however, since it occurs at the region 
where the curved shell joins the flat plate portion of the 
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hull (the boundary of the shell section). This particular 
region requires a separate analysis regardless of shell 
theory employed (bending or membrane theory) . In the 
former case, rather undefineable boundary conditions exist 
along with compatibility obstacles and in the latter case 
membrane theory never does acco^ant for conditions at the 
boxmdaries. Therefore, it is of no consequence that the 
radius of curvature expression given by the cycloid is 
invalid at the boundary of the shell since the shell theories 
are also in question at that region. A pertinent subject 
for investigation would be the study of the discontinuity 
stresses present at this junctioni Results could be 
correlated with this thesis which purports to treat only 
the bottom section of the ship not at or adjacent to the 
boundaries . 
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B. ORTHOTROPY OP SHIPS HULL 



As a ships hull is fabricated, the elements of 
cpnstruction (shell plating, longitudinal stringers, 
transverse frames) are arranged in such a manner that the 
resulting structural elastic properties ( extensional, bend- 
ing, shear and twisting rigidities) are not equal in all 
directions as they are for an isotropic structure. The 
varying of elastic properties with direction, or the aniso- 
tropy of the structure, is somewhat simplified by the 
symmetry of construction in the longitudinal and circumferen- 
tial directions. This special case of anisotropy where the 
elastic properties are different in mutually perpendicular 
directions is termed orthotropy and such structures are called 

r' 

orthotropic. It should be noted at this point that the basic 
material of construction is isotropic. The disposition of 
this material creates the anisotropy and therefore the 
structure can quite correctly be called constructionally 

orthotropic . 

The reaction of a structure to any given loading will 
be a function of its elastic properties. It would seem 
then that a realistic analysis could be made if the structures 
elastic properties and their change with direction are recog- 
nized. This can be accomplished by considering an ’’equiva- 
lent” structure without stringers and frames but whose elastic 
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properties and general behavior imder load resembles 
stiff iclently the original structure. 



The advantage of this approach is that the continuity 
of the structures behavior is retained and the resulting 
analysis is superior to one in which conservative attempts 
are made to account for the complex interactions existing 
between structural elements such as frames, longitudinals 
ahd shell plating. A critical question arises at this 
point - what is the characteristic mode of deformation of 
the stiffened hull? Does it deform under load to form a 
smooth wave length between bulkheads as in Figure 

or is all the deformation confined to the shell surface 
between stiffening members as in Figure 2l4.b ? 















— c 




% 




1 












'2-4 






24 b. 



FI go 2h. Stiffened Hull Modes of Deformation 
The answer is clearly a function of the number and size of 



the transverse frames, Lundgren (7) is satisfied to accept 
the former version of deformation (Figure 2li.a) j whenever the 
number of transverse frames between bulkheads equals or 

exceeds four. There are panel advocates^ however, who 

>• 

insist on the latter mode (Figure 24.b) as being the predomi- 
nant one which accounts for most of the strain energy of de- 
formation regardless of the number of frames. 

For the destroyer type having approximately i|- frames 
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between bulkheads (in the middle half of the ship), the 
author feels that the normal mode of deformation is probably 
a combination of smooth half wave length deformation with 
smaller panel deformation superposed on it as in Plgxire25 • 





1 1 























Fig, 25 Combined Modes of Stiffened Hull Deformation 
Perhaps many of those who can accept only panel deformation 

are misled by the physical appearance of ships hulls which 
have the so called "hxmgry horse” look or a permanent deforma- 
tion on all panels making the shell appear as though it were 
shrunk fit onto the framing system. These effects could 
certainly have originated from shrinking of the shell around 
faying surfaces welded to the hull or possible a dynamic 
loading effect of the ocean waves. 

The deformation of Figure25 (which is the authors 
version of the situation) illustrates that continuity of the 
structure between bulkheads must be considered in analysis. 

The extent of this significant structural section also requires 
a full consideration of shell curvature. One handicap of 
the gross section approach which is also obvious from Figure 
is that the localized stresses created by the superposed 
panel deformation would not be immediately determined. A 
final check of the solution should incorporate a method of 
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accoxmting for these localized stresses. 

The orthotropic rigidity coefficients (equation lij.) have 
been obtained and defined by the superposition of the elastic 
laws of the isotropic shell and the mating gridwork frame, 

A better appreciation for these may be achieved by an examina- 
tion of these coefficients on an actual hull. Figure 26 
illustrates the magnitude and variation of the coefficients 
for the DD 931 class destroyer. The data for these plots 
is in Appendix D 
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C. BOUNDARY CONDITIONS FOR USE WITH 



PARTIAL DIFFERENTIAL EC^UATIONS OF SHELL 



The solution of the developed partial differential 
equations ) will eventually require the stipulation 
of appropriate boundary conditions. These conditions will 
consist of a rational approximation of (a) the degrees of 
fixity existing at the bulkheads and at the junction of 
the curved shell and flat ships side and, (b) the magnitude 
and distribution of the loading which exists at the afore- 
mentioned boundaries. 

Fixity : 

The degree of fixity which exists at the boundaries 
of the shell could indeed be the subject for a separate 
and absorbing analysis. Vedeler (1) suggests a method 
for accomplishing this. It would be of questionable value 
however, to pursue this subject in any depth in this thesis 
since the immediate problem of determing a method of 
solution to the differential equation is the major obstacle. 
For a first approximation then,-, -bhe following conditions 
are postulated. At the junction of the shell and the bulk- 
head, the fixity can be considered to be 100^ clamped. The 
presuratlon is a natural result when one considers that the 
structure is a shell continuous over supports (bulkheads) 
with continuous and equal loading on each side of these 
supports. An examination of the bulkhead N^ stress 
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distribution in Figure 11 would also lead one to assTome 
that a full clamping condition prevails. The fixity 
condition along the longitudinal boundary is a bit more 
elusive. Once more, as a first approximation, it will be 
assxjmed as being simply supported and capable of movement 
in the athwartship direction. The condition can be 
represented by Figure 27 • 







1 / 

/ 

/ 


/ 




/ 


/. 




/ 


/ 




/ 


/ 










^ 








Figure 27 Degrees of Fixity at Shell Bovindaries 



Edge Loading ; 

To begin the estimate of edge loading and its dis- 
tribution, consider the bottom section and the total 

. OhI 

vertical force^due to the water (q) and the internal 
loading (w) . 







= 2 /( y ho COS0 - Yz cos$( - w)2H COS0 Ld^^ 



-TA 



= i|-HL y I (ho cos^p ■+ z cos^0 - w cos0) 



= ijUL 



Vt = HL 



fho cos^0 +H/2 (cos 20 4- cos^0cos20; - w cos^ d^ 

. L y J 

IT hp ^ 3T H - w •] 



y 



J 
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An alternate expression for could have been, 

\ - Wn - wt) dS^ 

where are normal and tangential components of the 

internal loading w. dS^is an increment of surface area 
projection on the horizontal plane. 

Substitution of Wj^ -= w cosO 

^ w sinO 

dS^ •= L cos^ d(6 

and Integration v/ill give same results obtained previously 
for V^. 

Having now obtained the total upward force, a 
determination must now be made of the distribution of the 
reaction forces which support V-^; and act on the shell 
boundaries. This is a statically indeterminate problem. 

In order to circumvent this obstacle and avoid an auxiliary 
analysis ^whose accuracy would be greater than justified;, 
the following approach is taken. 








Pig, 29 Analogy of Flat Plate Prismatic 
Loading to Hull Loading 
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Timoshenko (10 describes the behavior of a simply 
supported rectangular flat plate vinder a triangular 
prism distribution of loading. Figure 29. The portion' 
of this reference which is of interest is the nature of 
the reactive forces along the simply supported edges. 
Indeed it could be said that the vertical component of 
loading on the ships hull is a triangular prism distri- 
bution (with max value at the keelj if the hull were 
developed out into a flate plate. Figure 29* Because 
the ships loading is not on a flat plate but on a cylin- 
drical hull, no further resemblance or analogy will be 
made other than a consideratipn of the relative distribu- 
tion of the reactive forces and their relationship to the 
aspect ratio of the simple supports, Timoshenko tabulates 
^y max » ^x max various aspect ratio of simple support 
boiindaries. For the ship types being considered in this 
analysis, the pertinent ratio of B/L is approximately 1,3. 

2,73 X aspect ratio 



or 

2.73 L/B 



This value gives 

^y max 
^x max 
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Therefore, 



Vb . lb 
ft 



2LV^ + 2 BVg = 



2LV^ + 2 BVg 
and. 




HL V ("ir ho i- 3A'TT H - ^ ^ 



2.7 L/B 



Solving for the reaction forces gives. 




1.35 bhl^ Y t -rr ho 4 H - ij. 
2.7 4 -b 2 i|. 




0.5 BHL / ( iT ho ^ H - 
2.7 l2+b 2^ 1^. y' 




Fig, 30 Edge Loading on Shell Boimdarles 
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D. RIGIDITY COEFFICIENTS DATA 



The orthotropic rigidity coefficients of equation (II4.) 
can be written in the following form for 0.3* 

Extensional coefficients 

r ■ ^ 

E b£ 

D-v 0,33t 
E ■ 

^ __ 0*38t 

E 

Flexural coefficients 

^ _ ,09lSt2 ^ Ax°x^ 

E ^2 

_ .09lSt3 ^ 1|} 

E 

,06I|-lt3 ^ ,38 Jx 

E ’ '•^2 
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^0X - 

E 



- .0275t^ 

E 



Risldlty moments 



I 



I 



TABLE 3 DP 931 RIGIDITY COEFFICIENTS DATA 





0° 


10° 


30° 


45° 


72° 


0 

0 

0^ 

! 


(in) t 


.75 


.75 


• 625 


.4375 


.3125 


.375 


(in^) A^ 


19.8 


2.53 


2.53 


2.53 


2.04 


2.04 


(in) 


384 


38.il- 


38.4 


38.4 


28.8 


28.8 


(in) 


17.03 


6.73 


6.66 


6.57 


4.34 


4.37 


(in^) 


noil-. 


26.32 


26.32 


26.32 


7.61 


7.61 


(in^) Jx 


1215 


27.37 


27.37 


27.37 


8.81 


8.81 


(in^) A0 


13.5 


12.75 


8.7 


8.7 


5.61 


5.24 


(in) C0 


17.83 


16.38 


12 .66 


12.57 


9.66 


9.40 


(in^) l0 


859.5 


65o .5 


328.5 


328.5 


114.9 


107.8 


(ln> 4 -) 


931.5 


722.5 


344.5 


344.5 


121.7 


111.8 


(in) = 8I4. ~ const. 
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TABLE 14- DP 931 CALCULATED EXTJilMSlQML HIGIDITY 



COEPPICIENTS AT BULKHEAD 





0 ^ 


0 

10 


0 

30 


45“ 


0 

60 


0 

72 


80“ 




fl 3936 


+13341 


+ 9150 


+ 3200 


-3730 


44800 


148000 


tsi 

1 


•35805 


-34360 


-24950 


-I 5 y 380 ' 


-7360 


- 522^ 


- ni? 




0 


10753 


27910 


33920 


33800 


32690 


38850 


Ti 


16600 


16220 


13450 


10470 


7300 


5310 


4480 




1.215 


1.08 


0.58 


0.36 


0.308 


0.293 


0.88 


E E 


1.339 


1.19 


.638 


.396 


.339 


.322 


.969 


E 


0.4’ 


0.356 


0.191 


0.119 


0.102 


.096 


0.29 


^ A 

x0 

E 


0.461 


0.41 


0.22 


0.137 


0.117 


0.111 


0.334 
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table S dp 931 RIGIDITY COEFFICIENTS (AcUo-L) 





0° 


10° 


30 ° 


450 


720 


■90O 


s 


1.341 


0.891 


0.754 


0.547 


0.415 


0 .484 


“0 

ET 


0 .986 


0.977 


0.792 


0.585 


0.411 


0.475 


E 


0.248 


0.248 


0.206 


0.144 


0.103 


0.124 


E 


0.285 


0.285 


0.238 


0.166 


0.119 


0.143 


E 


178.8 


3.70 


3.64 


3.54 


1.60 


1.62 


H 

E 


61.45 


48.39 


20.63 


20.36 


7.56 


6.7 


Ks> 


0.0115 


0.0115 


0.0067 


0.0023 


0.0008 


0.0015 


E 


12.03 


0.298 


0.287 


0,276 


0.118 


0,120 


E 


4.237 


3.297 


1.576 


1.565 


0.552 


0.539 


Sx 

£ 


8.8 


0.44 


0.44 


0.434 


0.308 


0.310 


E 


2.87 


2.48 


1.32 


1.31 


0.64 


0.58 , 
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TABLE 6 DE 1033 CALCULATED BXTEMSIOML RIGIDITY 
COEFFICIENTS AT BULKHEAD 



0 


0“ 


10° 


30° 


45“^ ■ 


60° 


0 

72 


80“ 


- Nx 
ft) 


-»7725 




+5201 


+ 2180 


-750 


<-8000 


+110,000 




-2ij.,580 


-23650 


-17330 


-10880 


-5375 


-2463 


-1162 




0 


6463 


16910 


20,770 


21,250 


21,600 


27,200 


(tsi) fl 


15000 


14680 


12150 


9450 


6600 


4800 


4050 




0.538 


0.495 


0.321 


0.240 


0.198 


0.202 


0.635 


E E 


0.592 


0.545 


0.354 


0.264 


0.218 


0.222 


0.7 


Q 1 


0.1775 


0.103 


0.106 


0.0792 


0.0654 


0.666 


.21 


E 


0.20ii. 


0.188 


0.122 


0.091 


0.0753 


.0768 


0.242 
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TABLE 7 DE 1033 RIGIDITY COEFFICIENTS 



0 


o 

0 


10° 


& 

30 


0 

45 


0 

60 


0 

72 


0 

80 


t 


0.5 


0.5 


0.11-375 


0.4375 


0.375 


0.375 


0.3125 






2.81 


2.49 


2.20 


2.20 


2.20 


2.20 


b2 


26.9 


26.3 


27.5 


25.6 


26.3 


28.1 


26.2 


Ajf 


8.5 


7.3 


4.63 


4.57 


4.27 


4.03 


3.91 


E 


l.i|7 


0.657 


0.571 


0.566 


0.497 


0.489 


0.428 


E 


0 .61^4 


0.631 


0.531 


0.53 


0.46 


0.458 


0.387 


E 


0.165 


0.165 


0.144 


0.144 


0.124 


0.124 


0.103 


E 


0.19 


0.19 


0.166 


0.166 


0.142 


0.142 


0.119 


r ' \ 

ii 1 

b, = 90 = const, 1 ! 

1 1 

1 1 


1 

f 

f 


\ 

i 

t 
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